Light-induced breaking of symmetry in photonic crystal waveguides with nonhnear 
defects as a key for aU-optical switching circuits 



O 
(N 
> 

Q 

o 

(N 



O 



Evgeny Bulgakov, Almas Sadreev, and Konstantin N. Pichugin 
Kirensky Institute of Physics, 660036, Krasnoyarsk, Russia 
(Dated: October 27, 2011) 

We consider light transmission in 2D photonic crystal waveguide coupled with two identical non- 
linear defects positioned symmetrically aside the waveguide. We show that with growth of injected 
light power there is a breaking of symmetry by two ways. In the first way the symmetry is broken 
because of different light intensities at the defects. In the second way the intensities at the defects 
are equaled but phases of complex amplitudes are different. That results in a vortical power flow 
between the defects similar to the DC Josephson effect if the input power over the waveguide is 
applied and the defects are coupled. As application of these phenomena we consider the symmetry 
breaking for the light transmission in a T-shaped photonic waveguide with two nonlinear defects. 
We demonstrate as this phenomenon can be explored for all-optical switching of light transmission 
from the left output waveguide to the right one by application of input pulses. Finally we consider 
the symmetry breaking in the waveguide coupled with single defect presented however by two dipole 
modes. 

PACS numbers: 



I. INTRODUCTION 



p ^ , Symmetry breaking in a nonlinear quantum system is a fundamental effect caused by the interplay of nonlinearity 
• ■ with linear potential which defines the symmetry. It is commonly known that the ground state in one-dimcnsional lin- 
. 3 ', quantum mechanics is nodeless and follows the symmetry of the potential. However the self-attractive nonlinearity 

■ in the nonlinear Schrodinger equation breaks the symmetry of the ground state, replacing it by a new asymmetric 
I I state minimizing the systems energy. For example, the nonlinear Schrodinger equation in double- well potential reveals 

anti-symmetric ground state with variation of normalization of the state The phenomenon of the spontaneous 
' symmetry breaking in analog with the double-well potential are realized in a nonlinear dual-core directional fiber 0- 
^ . 1^- Spontaneous symmetry breaking was demonstrated recently by Brazhnyi and Malomed in a linear discrete chain 
, (Schrodinger lattice) with two nonlinear sites Q . They have shown as analytically as well as numerically the existence 
^ ^ of symmetric, anti-symmetric, and non-symmetric eigen-modes with eigen-frequencies below the propagation band of 

■ the chain, and that a variation of the population of modes can give rise to a bifurcation form one to another mode. 
The system has direct relation to photonic crystal (PhC) waveguides with two in-channel nonlinear cavities where the 
population of the cavities might be governed by external source of the light. 

Indeed, the phenomenon of symmetry breaking is studied in the nonlinear optics with injection of input power with 
the establishment of one or more asymmetric states which no longer preserve the symmetry properties of the original 
state (6l-[To| . In particular Maes et al. [ll|, [l^ considered the symmetry breaking for the nonlinear cavities aligned 
along the waveguide, that is a Fabry-Perot architecture close to the system considered in Ref. [8|. That system 
is symmetric relative to the inversion of the transport axis if equal power is injected on both sides of the coupled 
' cavities. The symmetry breaking was found also for the case of many coupled nonlinear optical cavities in ring-like 
architecture [isl [l^ . In the section II we write the equations of motion for the nonlinear optical cavities coupled with 
PhC waveguides by using an analogy of the two-dimensional PhC with quantum mechanics [15i] . As for an application 
we consider three simple PhC systems which undergo the symmetry breaking phenomena for variation of the light 
frequency or the input power. The first simplest system is two identical nonlinear defects positioned symmetrically 
aside the straight forward linear waveguide (section III). Each defect is presented by single monopole mode. We 
show two types of the symmetry breaking [l^, . In the first type the symmetry is broken because of different light 
intensities at the defects. In the second type of the symmetry breaking the intensities at the cavities are equaled but 
phases of complex amplitudes are different. That results in a vortical power flow between the defects similar to the 
DC Josephson effect if the input power over the waveguide is applied, and the defects are coupled. 

In section IV we consider as the phenomenon of the symmetry breaking can be explored for so called all-optical 
switching p^| - |2l| by use of the T-shaped photonic waveguide with two identical nonlinear cavities positioned sym- 
metrically. That system combines two systems. The first one is the Fabry-Perot interferometer (FPI) consisting of 
two nonlinear off-channel cavities aligned along the straightforward waveguide considered in Refs. [HI [H, H^l ■ As was 
shown in Ref. [11] there is a discrete set of the a self-induced bound (localized) states in continuum (BSC) which are 
the standing waves between off-channel cavities. In the second system two nonlinear cavities are aligned perpendicular 
to the input waveguide. As was said above there is the anti-bonding bound state in continuum (BSC). Here we show 
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that both types of the bound states might be important for the breaking of symmetry. All these phenomena agree 
well with computations based on an expansion of the electromagnetic field into optimally adapted photonic Wannier 
functions in two-dimensional PhC [l^, [l^l • 

Finally, in section V we consider the single nonlinear defect with two dipole eigen-modes which belong the prop- 
agation band of the PhC straightforward waveguide (section III). We demonstrate the symmetry breaking provided 
that the system is excited with equal powers from both sides similar to that Maes et al. has shown in the system of 
two coupled nonlinear cavities [llf, [l3 • 



II. BASIC EQUATIONS 

The light propagation in linear PhC is described by the Maxwell equations 

D{r,t) = eoir)E{r,t)- 

We take the light velocity to be equal to unit. However if there are defects with instantaneous Kerr nonlincarity, the dis- 
placement electric vector interior to the defects has a nonlinear contribution D{r, t) — eo{r)E{'r, t)+x^'^^ [E^i^i t)]'^E{'r, t) 
[H [li]. A substitution of the electric field in the form [E{f,t) = i[£;(r)e'"* -I- i?*(f)e-*'^*] into Eq. (P) and neglect 
by highly oscillating terms such as e^*" allows us to write the Maxwell equations in the same form as Eq. ([T|) with 

e{r) = eo(rl + \x^'\uj)\E{f)\^ E{f) + \x^^\u^)E\r)E* {r). (2) 

In what follows we consider the 2D PhC with arrays of infinitely long dielectric rods as shown in Fig. [ija) in which 
the electric field is directed along the rods while the magnetic field is directed perpendicular to the rods [in the plane 
of Fig. [Ha)]. Then Eq. © simplifies as follows [H] 

- eo(r) + -J'\Lo)\E{r)\''E{f). (3) 



There is a remarkable analogy of electrodynamics in dielectric media with quantum mechanics [15l.l27l|. In particular, 
if the nonlinear contribution to the dielectric constant is small we can use the well-known methods of quantum 

mechanical perturbation theory. Let |^) = _i. be the electromagnetic state in the PhC. Then the Maxwell 



equations ([T]) can be written as the Schrodinger equation — H\%1)) indeed with the Hamiltonian |15l |271 128| 



-iVx 

Because of the perturbation of the dielectric constant ^ the Hamiltonian can be presented &s H = Hi} + V where 
and 

^e(f)) e(r) eqW 

Let us introduce (following, for example, Refs. (28l. [29j) the following inner product for the unperturbed system: 

[e^{r)E*E' + H*H']d^r. (7) 
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which obeys the following normalization and orthogonality conditions for the bound eigen-states of the unperturbed 
Hamiltonian -ffolV'm) = '^mlV'm) 



(V'nIV'n') = ^ J [eo{T^E*„E„, + H*H„,]d^r = J eo(f)^:K'dV= <5„„,. 
Then the matrix elements for the perturbation calculated by use of these eigen-states are 



(8) 



(9) 



One can see that the matrix © is not Hermitian as was noted in Ref . [28| . The origin is that the unperturbed states 
obey the inner product ([7]) with the dielectric constant eo{r) while the eigen-states of the full Hamiltonian Hq + V 
obey the inner product with a different dielectric constant e{r). Respectively, the Hamiltonian H is non- Hermitian 
with the inner product (|H]). 

In order to avoid this problem we must use the inner product which is not tied to a specific choice of the dielectric 
constant. One way, given in Ref. 15[, is by using only the magnetic field for the state. Another way is to absorb the 
dielectric constant in the scalar product by a new function as F — ^ e{f)E. Then the inner product becomes 



[F*F + H*H ](fr. 



(10) 



The value = \ /[e(r)|£?p -I- \H\'^]<Pr is proportional to the energy of EM field which is important for the 

derivation of the forthcoming coupled mode theory (CMT) equations. That technique changes the Maxwell equations 
as follows: 



V X 
1 

7^1 



F 



= -H 
■.V X H = F. 



(11) 



The Hamiltonian takes the following form 





Ho 



-iV X 







Vx 



V 



V 



-iV X 6 



eir) 





Now the eigen-states of the full Hamiltonian can be expanded over the eigen-states \m) 
turbed Hamiltonian Hq where 

F 

V X 



Vx 



V X H,n = -ioJm\/^o(/)F„ 



Err 



(12) 



of the unper- 



(13) 



Then we obtain from © 



{m\V\n) = / d'rco' {r)6 



E*^ir)E„if). 



(14) 



One can see that the full Hamiltonian is Hermitian now. 

If the nonlinear defect rods are thin enough, the dielectric constant ^ can be rewritten as follows 



e,(x) ^ (.o + -x(^'(^)|i?(x)nE^(^-^.- 



(15) 
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Here j enumerates the defects, 6 = 1 inside the defect rod and 9 = outside. As was shown for the simple square 
lattice 2D PhC from thin GaAs dielectric rods [IJl the resonance spectra in the PhC waveguide are located in a rather 
narrow frequency domain. Therefore, we neglect the frequency dependence in the nonlinear susceptibility x^^H^) in 
the following. Assuming that the nonlinear contribution in Eq. (|15p is small compared to eo we obtain for the matrix 
elements (fT4| per unit length of the defect rods 

{m\V\n)^ L_^(3)(^^+^„)^ /■ d2^|i?(x)pi?,„(x)*£;„(x). (16) 

16eo j •''^j 

In order to find electric fields at the defects we must constitute a way to excite the defect modes. Here we consider 
that the EM field propagates from the left along the straight forward waveguide, interacts with the nonlinear defects, 
reflects back and transmits to the right. Then the transmission process can be described by the CMT stationary 
equations (30l - l33| 

[t^ - '^{LOmSmn + Vmn + «r„]A„j = iy/V^Ein. (17) 

n 

These CMT equations, in fact, are the Lippmann-Schwinger equation [s^, Issj 

[uj-Heff)-^ = iWE^n. (18) 

where the complex matrix H^f / equals 

Heff =H(i + V -iWW+, (19) 

the columns of the matrix W consists of coupling constants of the m-th eigen-mode with the p-th injecting wave 
^jTrap, and the column 5* consists of the mode amplitudes Am- The solution 4' is given by inverse of the matrix 

Lo — Heff where the matrix elements of the effective Hamiltonian Heff in turn depend on the mode amplitudes Am- 
In order to write the equations of self-consistency for the amplitudes at the defects we expand the electric held at the 
j-th defect over eigen- modes E{xj) = J2m ^mV'm(xj). That deflnes the equations of self-consistency after substitution 
into Eq. (fTT]). 

Finally, we present the transmission amplitude in the framework of the CMT [sil, H^l 

t = E^n~ W+-^. (20) 



III. LINEAR OPTICAL WAVEGUIDE COUPLED WITH TWO NONLINEAR OFF-CHANNEL 

CAVITIES ALIGNED SYMMETRICALLY 

Two identical nonlinear defects positioned symmetrically relative to the single straight forward waveguide is one of 
the simplest systems in which the breaking of symmetry occurs [l^ [l3| . The system can easily be realized in 2D PhC 
as shown in Fig. [T] (a). The system is symmetric relative to the inversion of the y axis, as shown in Fig. [Hb), and 
thereby supp lements the system in which two nonlinear cavities are aligned along the waveguide considered by Maes 
et al [lllll^l- That system is symmetric relative to the inversion of the x axis if equal power is injected on both sides 
of the waveguide. 

Let each defect supports a localized non degenerate monopole solution for the TM mode only, which has the electric 
field component parallel to the infinitely long rods (ill. [23|. Other solutions, (dipole, quadrupole, etc.) are assumed 
to be extended in the photonic crystal for the appropriate cavity radius and the dielectric constant j23 . [s^ and are 
thereby excluded from the consideration. Therefore, we have a two-level description for Hq with the eigen-frequencies 

LOs^a = UJQ±U (21) 

where u is the coupling constant u. We denote the corresponding even (bonding) and odd (anti-bonding) eigen-modcs 
as V's,a(x). Both modes for specific PhC are shown in Fig. [51 We pay attention that the frequency of the bonding 
(nodeless) mode is higher than the frequency of the anti-bonding mode with one nodal line. 

Next, we assume that the EM wave which propagates along the waveguide obeys the symmetry of the total system. 
Therefore the wave might be only symmetrical relative to ?/ — )■ —y or anti-symmetrical. Respectively, the symmetric 
wave could excite only the bonding mode while the anti-bonding mode would remain as a hidden mode. Other 
words, the anti-bonding mode is the bound state in symmetrical continuum ^33]. The same refers to the case of the 
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(a) 
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FIG. 1: (a) Two defect rods made from a Kerr media marked by filled circles are inserted into the square lattice PhC of 
dielectric rods with the lattice constant a = 0.5/im, the cylindrical dielectric rods have radius 0.18a and dielectric constant 
£0 ~ 11.56. The ID waveguide is formed by substitution of linear chain of rods by the rods with dielectric constant ew + eo 
marked by stars, (b) Schematic system consisting of a waveguide aside coupled to two single-mode cavities. The cavities are 
coupled each other via u. 
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FIG. 2: (a) Bonding (even) mode and (b) anti-bonding (odd) mode of two overlapped linear defects in the 2D PhC. The defects 
have the same radius as the radius of rest rods but different dielectric constant eo — 3. The frequency of the isolated defect 
equals 0.3593 in terms of a value 2nc/a. For the case of two overlapped defects shown here the frequency is split to be equal 
0.3603 (bonding) and 0.3584 (anti-bonding). 



anti-symmetric wave and the defect bonding mode. Therefore, for the linear case one can see the only resonance dip 
at Us, if the symmetric wave propagates along the waveguide. However due to the nonhnearity the light transmission 
acquires much more rich behavior because of spontaneous breaking of symmetry. 

Substituting two eigen-functions 'ips,a{'^) into Eq. and considering a radius of the defect rods are very thin 

compared to the characteristic scale of wave function we obtain 

{m\V\n) « -lax^^\uj^+tOn) ^ |i;(x,)PVm(xj)*V^„(x,), (22) 



where a is the cross-section of the defects. Finally, we obtain from Eq. (|22|) 

^ ~ ^ a.O0.0a(/l - h) UJa<t>l{h + h) ) ^ ' 

where 0^ — '0s(xi)Vct = ijjs{'^2)\f^i4'a — V'a(xi)Vo' = — '0a(x2)v^, and X2 are the positions of the defects in 
the two-dimensional PhC, and Ij = |i?(xj)p,j = 1,2 are the intensities of the electric field at the nonlinear defects. 
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In order to find electric fields at the defects we must constitute a way to excite the defect modes. Here we consider 
that the EM field propagates from the left along the waveguide, interacts with the nonlinear defects, reflects back, 
and transmits to the right. Then the transmission process can be described by the CMT stationary equations [sol - ls^ 
for the bonding mode amplitude Ag and the anti-bonding amplitude Aa 

[UJ - UJs - Xu!s(t>l{h + h) + iT]As - Xujo4>s(f>aih - l2)Aa = i\/TEin, 

-Xua(j>s(pa{h - l2)As + [LJ-OJa- ^i^a<l>l{Il + hMa = 0, (24) 

where only the bonding mode is coupled with the waveguide because of the symmetry. The equivalent model is shown 
in Fig. [11(b). 

The amplitudes Ag and Aa are given by inverse of the matrix given in the right hand of Eq. ((24|) whose matrix 
elements in turn depend on the intensities /i , /2 ■ In order to write the equations of self-consistency for the intensities 
at the defects Ij,j = 1, 2 we expand the electric field E{x.) at the thin j-th defect over eigen-modes 0s(x) as E{Kj) — 
X]m=s a ^m4'm{^j)- The cxpausiou can be specified as follows 

El = E{xi) = (PsA, + (PaAa, E2 = E{X2) = - (f>aAa (25) 

where symmetry properties of the eigen modes (jjmi^) were taken into account. Respectively, 

h = \(t>sA,+(l,aAa\^,l2 = {(psAs - <paAa\^ (26) 



which defines the equations of self-consistency after substitution into Eq. (|24p . In general they are rather cumbersome. 
Let us, first, consider the more simple case of the isolated defects so that the overlapping u can be neglected. Then 
the values of the eigen- functions at the defects are equal cpg = <t>a- Even in that simplified case the solution of Eqs. 
([24l) has cardinal features different from the case of the single nonlinear defect considered in Refs. [37 43]. These 
features are the result of the mutual interference of wave fiows refiected by the nonlinear defects. If det{Hef / — i^) 7^ 
the amplitudes of the mode excitement for the transmission can be easily found from Eq. (1241) as follows 

. _ iVr£;„,[cj-cjo(i + 2A/)] 



{lu ~ a;o(l + 2A/))2 - wg^s ^ ^r(w - wo(l + 2A/)) ' 



{oj - a;o(l + 2A/))2 -a;2A2 -)-ir(w - wo(l + 2A/)) ' 



— ^, , , . , o\ T\\2 . ,2 A2 I .--n/, . , . /i i o \ rW ' ' / 



where the values / = (/i +/2)/2, A = A(/i —I2) in turn depend on the mode amplitudes according to (|26p . Substituting 
these solutions into Eq. (1261) we obtain the following nonlinear equations of self-consistency 



h = 



-u;o(l + 2A/2)] 



2 



[uj - ujoil + 2A/i)]2[a; - a;o(l + 2A/2)]2 +V^[uj- uja{l + 2X1)] 



r ^ rg2^[^~^o(l + 2A/i)]2 

^ [a;-cJo(l + 2A/i)]2[a;~a;o(l + 2A/2)]2+r2[a;-wo(l + 2A/)]2- 

The solution of these equations gives the steady state for the transmission in the waveguide coupled with two nonlinear 
defects. Finally, we write from Eq. (j20p equation for the transmission amplitude: 



t^E,r,-VfAs. (29) 

The odd amplitude Aa does not contribute into the transmission amplitude because of the symmetry. 

In the forthcoming CMT calculations we fix the parameters of the CMT model as follows: = l,r = 0.01, A = 
—0.01. We consider the case of isolated defects u = Q,4's = (pa = ^ and the case of coupled defects with u = 0.01, (ps = 
l,(pa = l-l- Rigorously speaking these values u and (ps,<pa correlate with each other. However, in our model case, we 
disregard this correlation. 



A. Symmetry preserving solution 

We start with the solution Ei — E2 that preserves the symmetry. In this case the incident wave excites only the 
symmetric even mode As 

i^/VEin 
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FIG. 3: Frequency behavior of the intensities at the isolated defects u — 0. (a) Ei„ — 0.01, (b) Ein = 0.05. Here and in 
the forthcoming figures dashed blue line shows the symmetry preserving solution. Solid and dash-dotted red lines show the 
symmetry breaking solutions which has different intensities at the defects Ji and l2- Gray thick solid line shows a new phase 
parity breaking solution at which det{ui — Hefj) ~ 0. 




FIG. 4: Frequency behavior of (a) the bonding amplitude \As\ and (b) the anti-bonding amplitude \Aa\ of the model shown in 
Fig. m (b) for u = 0, = 0.05. 



as follows from Eq. (I?f)) with the only resonance frequency ujq{1 + 2X1) and the width 2r. The self-consistency 
equation for the symmetry preserving solution I — Ii ^ I2 simplifies 

I[{uj-LOo{l + 2XI)f+r^]=rEl. (31) 

That coincides with the equation of self-consistency for the single ofF-channel nonlinear defect obtained in Ref. [s^ . 
The solution of this cubic nonlinear equation is shown in Fig. [3] by dashed blue lines. The frequency behavior of the 
intensities inherits the linear case, as shown in the inset. With growth of the input power the resonance frequency 
shifts to the left because of the nonlinear contribution 2X1 as seen from Eq. (pO)). 

The frequency behavior of mode excitations \As\, |Aa| is shown in Fig. |4]by blue dashed lines. As seen from Fig. 
111(a) As has a resonance peak. Respectively, the transmission T — \t\'^/Ef^ has a resonance dip at the frequency 
a;o(l + 2X1) = a;o(l -f 2XEf^/T) as shown in Fig. [SJa) by the dashed line. The last equality follows from Eq. (I5T]) . 



B. Symmetry breaking solution 

For the transmission through the nonlinear symmetric media the symmetry might be broken [ol-fl^. Numerical 
solution of Eq. p4|. indeed, reveals the solution with Ji > I2, i.e., the nonlinearity gives rise to a breaking of the 
symmetry below (above) the critical frequency for A < (A > 0). The symmetry breaking solution is shown in 
Fig. El by solid lines for Ii and dash-dotted lines for 12- There is also the solution that differs from the former in that 
El ^ E2. If the solutions are stable, a choice of the solution happens incidentally, as it does for a phase transition of 
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FIG. 6: The difference between tlie intensities at the defects for u — 0: (a) Ein — 0.01 and (b) Ein ~ 0.05. Only tire symmetry 
breaking solution is shown. 



the second order in cooperative systems. As shown in Fig. |5J a value Ii — I2 or the odd mode amplitude Aa, indeed, 
might serve as the order parameter that characterizes the symmetry breaking. 

It is surprising that there is the frequency at which the intensity at one of the nonlinear defects turns to zero as 
shown in Fig. [3l According to Eqs. ([28| that occurs at the frequency 

lomp=u:q{1 + &\EIJT). (32) 

At this frequency we have As = Aa = Ei/2 = Ein/V^ in accordance with Eqs. p7)) and ((25)) . By substituting this 
equality into Eq. (|20|) we immediately obtain that the frequency (|32|) defines the position of resonance dip for the 
symmetry breaking solution. As will be shown, that result of full extinction of one of the nonlinear defects is observed 
in the PhC system as well [Fig. [T^ b')]. 

In Figs, m^a) and|4] (b) we show the frequency dependence of the even and odd mode amplitudes j^^l and \Aa\ 
respectively for Ei^ = 0.05. One can see that, first, the incident wave begins to excite the odd mode below uJc for 
A < 0, and, second, \As\ and \Aa\ show the bistability. The even mode Ag displays a resonance peak (solid line) with 
the resonance width twice less than the resonance width of the peak for the symmetry preserving solution (dashed 
line) . Correspondingly, the transmission in Fig. [5] demonstrates a narrow dip for the symmetry breaking solution. In 
order to understand that phenomenon let us consider the resonance poles of the even and odd amplitudes given by 
zeros of the denominators in Eq. (P7|) 

ZU2 = ^0(1 + 2A/) - y ± ^Joj^A^ - ^. (33) 

For the solution with A = we had the only resonance pole with the resonance half width F. As Fig. [S] shows there is 
the frequency domain roughly between 0.98 and 0.99 where ujqA > T/2 and where the resonance half-width is twice 
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less than T according to formula ([33]). Therefore, in this frequency domain we can expect the resonance dip to be 
twice narrower compared to the symmetry preserving solution with A = 0. 

The lesser the width of resonance, the more unstable the resonance [l^. One can thereby see that the bistability of 
the symmetry breaking solution is more profound in comparison to the symmetry preserving solution. The resonance 
peak in \As\ for the symmetry breaking solution terminates at that frequency where the odd mode amplitude \Aa\ 
arises as seen from Fig. IHb). Close to this frequency the amplitude Aa has a square root behavior typical for the 
order parameter in phase transition of the second order. The dependence of Aa on the amplitude of the incident wave 
demonstrates the same behavior [see below Fig. [3(b)]. 



C. Phase parity breaking solution 

At last, there is the solution that has equal intensities at the defects but nevertheless a symmetry is broken because 
of phases of the complex amplitudes Ei and £2- This solution refers to the special case of Eq. ([24]) when the 
determinant of the matrix to — H^ff equals zero, (i.e., the inverse of matrix does not exist). It occurs at 

1,^12^1, LU^Uail + 2X4>ll). (34) 

Then the solution of Eq. ([M|) for the even mode amplitude Ag is 

^ = (35) 

while Aa is undetermined yet. 

Let us take for a while, the defects to be linear. Then the second equation in (|34l) shrinks to the isolated point 
Lo = u>a- As given by the CMT equations ([24]) and as seen from Fig. [1] this odd mode has zero overlapping with the 
waveguide and Eq. ([M)) thereby defines the bound state in continuum (BSC) 34, 44-48]. The solution of the Eq. 
As\ 

1 with As given by Eq. ([55]) and arbitrary Aa is therefore a superposition of the transport solution and 
the BSC. 

For the nonlinear defects the situation changes dramatically. First, there is the whole frequency region w > Wq for 
A > or a; < Wa for A < where det{uj — Heff) = as seen from Eq. p4|) . Equation ([34)) thereby defines the BSC 
with eigen frequency in whole region as dependent on the BSC intensity. Second, the BSC can not be independently 
superposed to the transport solution for the nonlinear case. The BSC begins to couple with the incident wave and 
can not be defined as the bound state if Ein 7^ 0. 

(i) Let the defects be isolated; (i.e. u = 0,(f>s — (jya — ^- On the one hand, we obtain from Eq. (|35|) 



As = El±El ^ E,.J^, (36) 



according to Eq. (j25p . That is the bonding mode amplitude is constant over the frequency as shown in Fig. [D (a) by 
the gray thick solid line. On the other hand, Eq. (IMl) directly shows that the intensities at the defects do not depend 
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on Ein, 

Since \Ei \ = \E2\ ~ VT the only way to satisfy Eqs. ([36]) and ([37| is to consider that the amphtudes at the defects 
are Ei = ^/l cxp{i6),E2 = \fl exp(-i6'). That is illustrated in Fig. [S^a). With the use of Eqs. dM]) and dSZ]), we 
obtain 

cos^e^-?^^. (38) 

For Ein — ^ we have the following limits: Q — > 7r/2, E^ — > iVT, E2 -i\/T, + i?2 as seen from Eq. ([55]) . As 
soon as Ein 7^ the defects amplitudes are seized to oscillate in fully anti-symmetric way as shown in Fig. [HJa). We 
emphasize that phase difference W has nontrivial behavior if the defects are nonlinear (A 7^ 0) and the incident wave 
is applied {E^ 7^ 0) as follows from Eq. (|38|) . For the symmetry preserving solution = (dashed line in Fig. |9]), for 
the symmetry breaking solution 20 = or tt (solid line in Fig. ^ while for the present solution the phase difference 
20 behaves as an order parameter (gray thick dashed line in Fig. |9l) similar to Aa shown in Figs. SJ^b) orlTJb). 

We define the present solution of the CMT equations (fM)) with the zero determinant detiuj — Heff) = as the 
phase parity breaking solution . It exists for u! < luq + 2XEf^/T for A < 0. Knowledge of the phase 9 allows us now 
to find the anti-bonding amplitude 

Aa^ (El- E2)/2^iVl sine. (39) 

The frequency behavior of the even and odd amplitudes are shown in Figs. Hl^b) andjTJb). Finally, by 

substituting Eq. I\36\i into Eq. (|29| we obtain t = for the phase parity breaking solution as shown in Fig. [Sfa) by 
gray thick dashed line. 

(ii) Coupled defects. For the PhC structure shown in Fig. [IJa) the coupling between the defects u is rather small 
compared to the coupling between the waveguide and defects Vt. Nevertheless, an account of the coupling between 
the defects has a principal importance as will be seen below. As was given earlier, the parameters of the coupled 
defects are specified as follows u = 0.01, (j)s = 1, 0a = l-l- 

A substitution of Eq. (jM]) into Eq. ([35]) gives 

^ _ UJQ+U VVEin _ AsO ^^p^ 

Wo(l — a) + u(l + a) uj — ujr + iTr UJ — cj,. + iTr 

where 

(1 - a)WsWa 



(41) 

= r— ° , (42) 



a;o(l 




- a) ' 










-a) + u{l H 


- a) ' 




V^EinUJa 





ujq{1 - a) + u{l + a) ^ ' 
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FIG. 9: Difference between phases of the amplitudes Ei and E2 for u — 0.01 as a function of (a) the frequency for = 0.05 
and (b) the amplitude of the incident wave for oj = 0.95. Dashed blue line shows the symmetry preserving solution, solid red 
line shows the symmetry breaking solution, and gray lines show the phase parity breaking solution, u — dashed and u = 0.01 
solid. The BSC point is shown by open bold circle. 



a — 4>1/4>1,- Therefore for the coupled defects the amplitude Ag acquires typical Bright- Wigner resonance behavior 
in which the nonlinearity is excluded. Respectively, a substitution of the solution PO)) into Eq. (|29l) immediately 
results in the transmission having the resonance dip at the frequency ujr with the half width which depends on 
ratio a and u. That result is shown in Fig. [S^b) by the gray thick line. If u — > 0, 0a — ^ </'s, — ^ 1 the frequency of the 
resonance dip goes away, and F^ — >■ 00; that is, the resonance at the phase parity breaking solution disappears, and 
the corresponding transmission tends to zero as seen from Fig. [SJa). 
Equation fixes intensity at the defects 

which is similar to the former case given by Eq. (P7)) . On the other hand, we have according to Eq. (pS)) Ei + E2 ~ 
As/ 241 s where As is given by Eq. (PH)) . A graphic illustration of the solution of this equation with modules of 
Ej,j = 1,2 fixed by Eq. dH]), is shown in Fig. Hljb). By presenting Ei = ^/l exp(i(/3 + 0)) and E2 ^ exp(i(^ - 9)) 
we obtain from Eqs. (|40| 



cos^ 9 



tun 13 =^^. (45) 



The behavior of the phase difference 29 on the frequency or the incident wave amplitude Ein for u ~ 0.01 is shown in 
Fig. El 

However, the most remarkable feature of the phase parity breaking solution for u 7^ is related in a current circulated 
between the defects. When the phase difference 29 exists between two quantum dots (QD) or superconductors, 
connected by a weak link, a tunneling or Josephson current J = Jq sin 29 will flow between them. The value of 
the current Jq is proportional to the coupling constant between QDs or superconductors 49]. In order to explicitly 
write the expression for a current flowing between defects we use the Green function approach developed in Refs. 
[13,131,113 for the 2D PhC of dielectric rods with the dielectric constant eo- The PhC holds the ID cavity (waveguide) 
and two OD defects (nonlinear cavity rods) as shown in Fig. [1] (a). Then the dielectric constant of full system e(x) is 
a sum of periodic perfect PC and cavity-induced terms e(x) = epc(x) -I- Se{x\E), where Se{x\E) = eH'(x) -I- ed{x\E) 
is contributed by the waveguide and the two nonlinear defects: 

00 

ed{x\E) =ew ~ + E (46) 

n— — 00 j — 1,2 

Here 9 = 1 inside the cavity rod and ^ = outside, and the nonlinear contributions ej are given by Eq. (|15p. Then 
the TM electric fleld directed along the rods of the PhC E{yi,t) = £;(x)e*'^* is satisfled the integral equation 

i?(x) = d^yGix,y\cu)Seiy\E)E(y) (47) 
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where G(x, y|a;) is the Green function of the ideal 2D PC of the rods which was calculated in Ref. [Hoj] for the square 
lattice PhC. If the radius of the defects rods is sufficiently small in comparison to the wavelength of the EM wave, we 
can write Eq. (j47|) as the discrete nonlinear equation [33, [SO] 

En^Y. Jn-MSernE^ (48) 
m 

where Jn-m(w) = cr^G(Xn, x^\uj), a is the cross-section of the rods, and n, m runs over sites of the defects [marked 
by stars and filled circles in Fig. [T] (a)]. 

We use the nearest-neighbor approximation and write P5|) as the tight-binding linear chain coupled with two 
nonlinear defects 

- Muj)]En = Jl(-B„+1 + S„_i) + Sn,07^{S€iEi + Se2E2), 

[1 - SeiJoiuj)]Ei = J2iwEo + J4(5e2-B2, 

[1 - 5e2M^^)]E2 = J2ewEo + JiSeiEi. (49) 

The model is shown in Fig. [TU] and consists of a linear infinitely long tight-binding chain presented by amplitudes 

En whose spectrum is given by dispersion equation Jo{uj) = 2Ji cos A:, and two nonlinear defects presented by 

amplitudes 0i, (j)2- The coupling J2 connects the defects and tke chain and the coupling J4 connects the defects. 



— o — o — o 

J1 J1 J1 



FIG. 10: Tight- binding version of the PhC system shown in Fig. [T](a): J2 couples the chain and the defects and J4 connects 
the defects to each other. 

By multiplying Eq. (|49| by Eq = t* and subtracting the complex conjugated terms one can obtain the value of the 
power flow current flowing between the chain at the "0"-th site and the defects enumerated as j = 1, 2 as follows 

ja^i,2 ^ iwJ2lm{tEl 2)- (50) 

Similar manipulations with the cavity's amplitudes give the current between the defects 

ji^2 = Ji6eIm{EiE;) = J^Sel sm(20). (51) 

It follows also that the current from the "-l"-th site to the "0"-th site of the chain coincide with the current from 
the "0"-th site to the "l"-th one. Therefore the currents (ISUl) and ([5T|l coincide also in accordance to the Kirchhoff 
rule. Thus, the input power induces vortical current between the waveguide and defects via the couplings J2 and 
J4. The current is excited by the incident wave provided the defects are nonlinear. Thus, our analysis shows that 
the symmetry can be broken not only because of different intensities at the defects but also by a circulating current 
between the defects although the intensities at the defects are equal. This model result of the Josephson like current 
between the defects with different phase is reflected in computations of the Poyinting vector in the PhC structure as 
will be shown below. 



D. Stability of solutions 

Furthermore, we studied stability of different solutions by standard methods given for example in Ref. ^2^. The 
stability of the solution can be found from the temporal CMT equations 

ihs = [ijJs + \ujs(j)l{Ii + h) - iT]as + \ujo(t)s(t)a{h - h)aa + VFE'^e"*'^*, 
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CD 

FIG. 11: Regions of stability of the solution. The symmetry preserving solution is stable everywhere except interior of the 
closed region shown by solid blue line. The stability of the symmetry breaking solution is shown by red, while the phase parity 
breaking solution is shown by gray. The parameters are u = 0.01, tJo = 1, T = 0.01, A = —0.01, (j>s = 1, 0a = 1.1. 

icia = Xu!Q(f>s(l)aiIl - l2)as + [l^a + Xu!a(f>l{Il + ^2)]aa- (52) 

By presenting as{t) — {As + ^s(^))e^*"*, aa{t) = {Aa + ^a(i))e~'"* with As, Aa as the steady state obeying the 
stationary CMT equations and |Cs(^)| ^ \^s\, \£,a{t) \ ^ \Aa\ we obtain the hnearized time-dependent equations 
for complex ^s,^a 



( Reiis) \ 




/ Reiis) \ 


Imiis) 




Imiis) 


Reiia) 


Reiia) 






V Imi^a) J 



(53) 



Their stablHty is determined by eigen values of the matrix L which is time independent. The results of our calculation 
of stability are presented in Fig. [S] which shows that the stability of the phase parity breaking solution appears if 
only the defects are coupled and 0s 7^ fpa- We collected the results of stability of all three solutions in Fig. [11] in the 
form of phase diagrams in plane of the incident wave amplitude and the frequency. One can see that the phase parity 
breaking solution is stable in some small area of the phase diagram. 



E. Numerical calculations in photonic crystal 

We numerically solve the Maxwell equations ([T]) for the TM mode in the PhC with defect nonlinear rods by expansion 
of electromagnetic field over maximally localized photonic Wannier functions [1^ [U, H^l • The square lattice 2D PhC 
has the same parameters as given earlier [see Fig. HJa)]. For the case of isolated linear defects with the same radius 
as the radius of host rods and the dielectric constant = 3 their eigen frequency ujq = 0.3593 in terms of 27rc/a. 
Overlapping of the defect's monopole modes gives rise to splitting of this frequency cOg ~ 0.3603, cJq = 0.3584 as 
numerical computation of equations ([T]) gives. Respectively we obtain that the value of coupling u = —0.001. The 
corresponding bonding and anti-bonding modes for the nearest distance a between defects were found in Ref. [5l| . 
For more distance 4a they are shown in Fig. [21 By the normalization condition ([S]) the heights of the amplitude 
modes at the defects equal (j)s = 0.5569, 4'a = 0.6179. Let us evaluate the dimensionless nonlinearity constant A. We 
take in numerical calculations the incident power per length of order lOOmW/a which corresponds to the incident 
intensity /q = lOOmW/a^. For chosen PhC lattice with period a — O.biim we obtain that the incident intensity equals 
Q.QAGW/cm^. With the use of e = eo + 2^60712/0 we obtain 

A = -27^712/0. (54) 

We take the linear and nonlinear refractive indexes of the defect rods to be, respectively, no — ^Jeq = \/3, "-2 — 
2 X 10~^^CTO^/VF. By substituting all of these estimates into (|5^ we obtain A ~ —0.9 x 10~^ which is close to that 
used in the CMT consideration. Finally, we estimate the coupling of the defect mode with the propagation mode of 
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CO 3/271 C CO 3/271 C 

FIG. 12: Self-consistent solution for (a) the intensities of the EM field at the nonlinear defects and (b) transmission spectra 
in the PhC structure in the PhC structure shown in Fig. [T] The parameters of the PhC and defects are given in Fig. [5] The 
input power per length equals lOOmW/a. n2 = 2 x 10"^^cm^/M^, A = -0.009. 




FIG. 13: Absolute value of the EM field solution for (a) the symmetry preserving solution ija/2nc = 0.355, (b) and (c) the 
symmetry breaking solution for ijja/2nc — 0.355 and ijja/2nc = 0.358 respectively. The EM wave incidents at the left of the 
waveguide. 

the PhC waveguide VT. There are many ways to calculate F using for example Refs. [l^ [H, [s^. In the present 
paper we estimated F numerically by using the following approach. We took the single hnear defect aside the PC 
waveguide as shown in Fig. [Tfa), and calculated the transmission spectra. By the resonance width of the spectra we 
evaluated F = 0.00185. 

The self-consistent solutions are presented in the form of the intensities in Fig. [T2l a). which are similar to the 
CMT results shown in Fig. (Ha). Also, one can see three solutions in the transmission shown in Fig. I12f b). as was 
found in the CMT model for the transmission shown in Fig. [5] Fig. [T3] shows the EM field (the absolute value of 
the electric field) for the symmetry preserving solution (a) and for the symmetry breaking solution (b) and (c). In 
the latter case one can see that the field is strongly different at bottom and top. Moreover Figs. [3l and fT2l show that 
there is a frequency at which the intensity of the EM field might be zero at the bottom defect. Indeed, Fig. [TST c) 
demonstrates this case. 

In agreement with the model consideration Fig. [HI shows that current flows (the Poyinting's vector patterns) are 
strongly different for the different solutions. For the symmetry preserving solution we have laminar current flow over 
the waveguide with excitation of two current vortices around each defect. The laminar flow over the waveguide and the 
vortical flows around defects are well separated. The whole current pattern is symmetrical relative to the symmetry 
transformation y — > —y. The picture has a similarity with ballistic electron transport in waveguide coupled to an 
ofi-channel quantum dot '54']. For the case of the symmetry breaking solution there is a current vortex inside the 
waveguide complemented by two vortices near each defect, as shown in Fig. [l3](b). The circulation in vortical flow 
around the upper defect is opposite to the circulation around the bottom defect. The vortical flow in the waveguide 
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FIG. 14: Current flows for the symmetry preserving solution which inherits hnear case (a), the symmetry breaking solution 
(b), and (c) the phase parity breaking one at au}/2-Kc = 0.35. Bold open circles mark the nonlinear defects. 



and the vortical flows around the defects are well separated for both solutions. In the third case (c) for the phase 
parity breaking solution one can see the current vortex in the waveguide and single vortices around the defects are 
mixed. Nevertheless because of the continuity equation in the space beyond of the nonlinear defects V j (a;, y) = 
the vortical flows around the defects and in the waveguide can not cross. 



IV. THE T-SHAPE WAVEGUIDE COUPLED WITH TWO NONLINEAR DEFECTS 

One of the most ambitious goals in nonlinear optics is the design of an all-optical computer that will overcome 
the operation speeds in conventional (electronic) computers. Vital in this respect is the design of basic components 
such as all-optical routing switches and logic gates. It is believed that future integrated photonic circuits for ultra 
fast all-optical signal processing require different types of nonlinear functional elements such as switches, memory and 
logic devices. Therefore, both physics and designs of such all-optical devices have attracted significant research efforts 
during the last two decades, and most of these studies utilize the concepts of optical switching and bistability. One 
of the simplest bistable optical devices which can be built-up in photonic integrated circuits is a single cavity coupled 
with optical waveguide or waveguides [15;, i55| . 

The concept of the all-optical switchingis based on a discontinuous transition between the symmetry breaking 
solutions by a small c hang e of the inp ut [19l | . Many of these devices employ a configuration of two parallel coupled 
nonlinear waveguides (l8l [20l [56l458j |. Recently Maes et al demonstrated the all-switching in the system of two 
nonlinear micro-cavities aligned along the waveguide TT] by the use of pulses of injected light. In the present section 
we use similar approach to demonstrate the all-switching effects in the T-shaped waveguide coupled with two nonlinear 
micro-cavities [sol, [60|. 

We consider the PhC shown in Fig. [15] with the same parameters as given in Section III: the lattice constant 
a — Q.bfjim, the cylindrical dielectric rods have radius 0.18a and dielectric constant e = 11.56. We substitute two 
defect rods of the same radius as shown in Fig. [Tnimade from an instantaneous Kerr media with the nonlinear refractive 
index n = uq + rt2^o where no = -v/S and ?i2 = 2 x 10~^^cm^/VI^. The corresponding equivalent configuration of the 
T-shaped waveguide with two nonlinear defects is presented in Fig. [12] 

We start with the position of the defect rods shown in Fig. [15] (a) which have strong coupling with the output 
waveguides 2 and 3, and negligibly weak coupling with the input waveguide 1. We consider a light given by the 
amplitude 5i+ is incident into the waveguide 1 and outputs into all three terminals as shown schematically in Fig. 
1161 The outgoing amplitudes are labelled as 5*1-, 5*2- and 5*3-. Each nonlinear optical cavity is assumed to be given 
by single mode amplitudes Aj, j = 1,2 and coupled with the guides 2, 3 via the coupling constant 7 shown in Fig. 
[TSIbv dotted lines and with the guide 1 via the coupling constant F. 

We consider that the defects are not coupled. Therefore the eigen frequencies of the system of the defects equal the 
monopole eigen frequencies of the isolated defects ujj = wq -I- A|Aj p, j = 1, 2 shifted because of the Kerr effect. Then 
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FIG. 15: T-shaped waveguide with two nonlinear defect rods. The cases (a)-(c) differ by the positions of the nonlinear defects. 





B "2-^' 



FIG. 16: GMT model of the T-shape photonic crystal waveguide coupled with two nonlinear defects shown by filled bold circles. 
The defects are coupled with input waveguide 1 via the constant F shown by dashed lines and with the output waveguides 2,3 
via the constant 7 shown by dotted line. Separately each connection is shown with corresponding light amplitudes. 



the equations fTT]) will take the following form 



(cj - cjo - A|Ajp + ^7 + iT)Ai + irA2 = i\/f(S'i+ + cri_e^^) + ^^7^2-6 
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(55) 



Here phases 9 and (j) as shown in Fig. [16] are the optical lengths through which light goes between the connections^ 
These equations are to be complemented by the equations for light amplitudes at each connection A, B, and D 



(w — cji + 17 + iT)Ai + iTA2 = i\/T{Si+ + cri-e'^) + i^a^^e^' 



(56) 



These CMT equations are to be complemented by the equations for light amplitudes at each connection A, B, and D 

ai+e-'« = 5i+-Vr(Ai+^2) 
Sx- = cTi-e"' - Vr(Ai + yl2) 
^2- -fT2-e"^- 77^1 
53- = as-e'^ - 77 A2. 
a2+e-*^ = -V7^i- 

(57) 
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TABLE I: Parameter sets of the CMT for PhC T-shaped waveguide shown in Fig. [15] 



Type of structure in Fig. 1151 


(a) 


(b) 


(c) 


r (in terms of a/2-Kc) 
7 (in terms of a/2-Kc) 
u)o (in terms of a/27rc) 


0.0002 
0.03093 
0.3609 


0.0002 
0.00189 
0.365 


0.00189 
0.00002 
0.3596 



The T-connection at the C point connects ingoing and outgoing ampUtudes by the S-matrix as foUows 

(Ti_ \ I a b b\ I (Ti+ \ 

CT2- \ = \ b c d \\ n2+ \. (58) 

(T3_ / \b d C J \(T3+ J 

In particular, the solution of the Maxwell equations for the T-shaped waveguide without defects gives the ma- 
trix elements of the S-matrix dMl) a = -0.3547 + 0.308i,6 = 0.6 + 0.173«,c = -0.4319 -I- 0.2271i,d = -0.568 + 
0.2225i at uja/2Trc — 0.35. Eqs. ([57| . and (|58p form a full system of equations for 11 amplitudes 

Ai, (71+, (Ti-, 0-2+, cr2- 1 fs+j cTa-i 'S'2-, 'S's-. Substituting S'i_|_ = £^j„e*"*,Ai 2 — j4i 2e~'"* we obtain after 
some algebra the following stationary CMT equations 

- Heff) ( ) = '^"'^ ( 1 ) ' ^^^^ 

where 

H^ff=[ -^H L02-^G)' (^^^ 



G = 7 + r(l + ae2*») + -fde^^'f' + y/^{b + c)e'^+"^, 
i7 = r + V7r(6 + c)e**+*'^ -I- Tae^*^ -f 7(ie2*"^)' 

i^ = \/r(l + ae2^'')-H776e^(^+'^). (61) 

We can calculate all parameters which are necessary in Eqs. (|59p . ((60|) . and (jFTj) . For the light transmission in 
straight waveguide coupled with the single linear off-channel defect positioned at different positions we able to extract 
the coupling constant of the cavity with PhC waveguide T and the eigen frequency of monopole mode wn . The results 
are collected in Table HI The limiting case of the T-shaped waveguide with F = is considered in Ref. [53| with results 
qualitatively close to those shown in Fig. 1191 Here we find the solutions with the substitution of concrete parameters 
listed in Table U and present the output light transmission to the left and to the right waveguides. We start with the 
case shown in Fig. [15] (a). 

Moreover we present real dispersion curve Lo{k) calculated for the straight PhC waveguide to find as the optical 
length given by the phase 9 ov (p depends on the frequency w. This curve is shown in Fig. [T7] by solid curve for 
the parameters of PhC given in Fig. [TJ In the vicinity of the BSC frequency Wco/Sttc = 0.3402 we approximate the 
dispersion curve as linear to obtain 

phase _ 17.68(^ - ujc)a , , 

^ ^ + r, ■ l"^J 

TT ZTTC 

parameters for the phases. 

The case of the T-shaped structure shown in Fig. [TSja) has an analogy with the Fabry-Perot interferometer (FPI) 
comprising two off-channel nonlinear cavities. As was shown in Ref. [22[ there is a series of the self-induced BSCs 
which are the standing waves between the off-channel cavities. It was shown that the BSCs exist for any distance 
between the off-channel defect rods because of their nonlinearity [22]. Similar BSC solutions are expected to exist 
in the present case of the T-shaped waveguide coupled with two off-channel cavities shown in Fig. ll5t aV However 
these solutions must be the standing waves with nodes at the point of connection of the waveguides. Therefore in 
the present case the BSCs are the anti-symmetric standing waves. Fig. [TS] shows one of these waves in the T-shaped 
waveguide with two nonlinear defect rods with the eigen frequency ujca/^irc = 0.3402 which satisfies the condition 

?>k{LOc)a = TT. 

For the linear case the BSC has zero coupling with the symmetric EM wave which inputs in the waveguide 1. 
However as was considered in section II nonlinearity gives rise to the important effect of the excitation of BSC by 
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wa/27tC 



FIG. 17: Frequency behavior of the optical length (phase) shown by solid line for the straight PhC waveguide which is fabricated 
by removing one row of dielectric rods. The parameters of the PhC are listed in figure caption of Fig. \T\ The BSC point is 
marked by star. 
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FIG. 18: Solutions of Maxwell equations for the case shown in Fig. [15] which demonstrate (a) the anti-symmetric FPI BSC 
(anti-symmetric standing wave between nonlinear off-channel defects marked by stars) with frequency u}c.a,/2'KC — 0.3402 for 
zero input power P = and (b) as this BSC mixes with incident EM wave to give rise to the symmetry breaking because of 
nonlinearity for P — 1.94VK/a. Yellow stars mark defect rods. 



the transmitted wave, (i.e., an interaction between the anti-symmetric BSC and the symmetric transmitted wave). 
As a result the total solution lacks the mirror symmetry of the PhC structure shown in Fig. [121 That is one of the 
scenarios of the symmetry breaking. Note, that Maes et al. [ll| have already reported the symmetry breaking in 
the FPI. In order for the FPI with two off-channel nonlinear cavities to have the mirror symmetry equal input power 
must be applied to both sides of the FPI [ll| . In our case of the T-shaped waveguide this symmetry is achieved by 
application of the input power via the additional waveguide positioned at the center of the FPI. 

Figure [TOT a) demonstrates the solution for the light intensities of the cavities Ij — c\Ej\'^/8Tr,j — 1,2 with broken 
mirror symmetry where Ej are the amplitudes of the electric field in thin defect rods for two values of input power. 
The solution converges to the BSC point marked by the star if the input power limits to zero. At this point the 
symmetry is restored. As a result for the input power P ^ the light transmission from the input waveguide 1 to 
the left waveguide 2 differs from the transmission and to the right waveguide 3 as seen from Fig. I19f b). Moreover 
the difference Tj^ — Tji crucially depends on the frequency in the vicinity of the BSC point. Figure ?? demonstrates 
the solution of the Maxwell equations for the z-component of electric field (the scattering wave function) breaks the 
mirror symmetry because of mixing the symmetric input wave with the anti-symmetric BSC as shown in Fig. I18f b). 
The intensity of the defect's modes is centered around the BSC intensity which is rather large as shown in Fig. [T9i;a) 
while the incident light intensity P ~ 1 • W/a. Thereby we have chosen exponential scaling for the solution of the 
Maxwell equations presented in Fig. I18f b) in order to distinguish waves in waveguides. 

As was shown in the framework of the CMT there might appear an additional branch of a loop shape for the 
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FIG. 19: (a) Frequency behavior of light intensities at the defect rods given in terms of W/a^. Dotted blue and dash-dotted 
red lines correspond to the input power equal to OASW/a, solid blue and dashed red lines do to 1.92W/a. Only those solution 
for intensities is shown which breaks the symmetry, (b) The frequency dependence of the transmissions to the left Tl (blue 
solid lines) and to the right Tr (red dashed lines) for the PhC T-shaped waveguide shown in Fig. [TJ] (a). The BSC point is 
shown by blue star. 
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FIG. 20: Frequency behavior of (a) the light intensities at the defect rods (in terms of Jo — W/a^) and (b) transmissions Tl 
and Tr for the PhC T-shaped waveguide shown in Fig. [TS](b). The input power equals 5.57VK/a. The symmetric solution is 
shown by solid gray line which inherits the linear case. The symmetry breaking solution because of the mixing the symmetric 
transport solution with the anti-symmetric BSC is shown by solid blue and red lines. The next symmetry breaking solution 
because of a bistability of the transmission in each output waveguides is shown by dotted brown and green lines. Zooming of 
loops are shown by arrows in (a). 

symmetry breaking solution with growth of the input power [s^ . The numerical results for the solution of the 
Maxwell equations completely agree with these model results as seen from Fig. [20la) by dotted brown and green 
lines. The loops shown by dotted lines in Fig. [501 are the result of individual instability that arises for transmission 
in the left or right waveguide coupled with the left or right nonlinear off-channel cavity ^59.] . Loops in the intensities 
are reflected in the loops in the transmission for the transmission to the left and for the transmission to the right as 
shown in Fig. I20r b) by dotted brown and green lines respectively. 

Bistability of the light transmission in the PhC waveguide coupled with nonlinear optical cavity crucially depends 
on the coupling: the smaller a coupling the less input power is needed for bistability [l^. For case (a) in Fig. [15] 
the coupling is rather large to observe bistability in the transmission. However case (b) has the sufficiently smaller 
coupling as one can see from Table |T1 As a result case in Fig. ITST b) gives rise to additional loops as shown in 
Fig. [50] for larger input power b.blW/a. Figure ^T\\ shows the wave function for the symmetry preserving solution 
(a) and the symmetry breaking solution (b) and (c). Cases (b) and (c) differ by the frequency. One can see from 
Fig- [2lT a) and Fig. [2lTb) that for the symmetric solution the transmission excites the cavities weakly, while for the 
symmetry breaking solution the defects are strongly excited because of mixing the injecting symmetrical wave with 
the anti-symmetric FPI BSC. Figure I^TT c) shows that for the frequency in the loop domain ujajl-Kc — 0.3442 the first 
nonlinear cavity is excited much more than the second one that is correlated with the outputs. 
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FIG. 21: The EM field solution in a scale real{Ez) exp(— jiJ^I) for (a) the symmetry preserving solution for ija/2nc = 0.3442, 
(b) the symmetry breaking solution caused by the BSC for uia/2vc = 0.3388, and (c) the for u)a/2-KC = 0.3442. Yellow stars 
mark defect rods. 
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FIG. 22: Frequency behavior of (a) the light intensities at the defect rods and (b) transmissions Tl and Tr for the PhC T- 
shaped waveguide shown in Fig. I15f c). The input power P — 0A8W/a. Only those solution is shown in (b) which demonstrates 
different outputs in the left (red dashed line) and right (blue solid line) waveguides. 

The T-shaped waveguide coupled with two nonlinear cavities shown in Fig. [15] is remarkable in that it allows the 
limit to the FPI [case (a)] with the FPI BSC in the form of the standing waves between two off-channel defects [lH as 
well as the limit to case (c) in Fig. [15] with the BSC in the form of the anti-bonding defect's state. Patterns of such 
anti-bonding BSC in PhC straight waveguide coupled with two cavities positioned perpendicular to the waveguide 
are shown in Fig. [5] 

The nonlinearity gives rise to mixing the anti-bonding BSC with the wave transmitted over waveguide 1. For the 
linear defect rods this state would be the perfect BSC. For the nonlinear cavities mixing this anti-bonding BSC with 
symmetric input light leads to the breaking of the mirror symmetry to give rise to the breaking of symmetry in the 
input waveguide. Then for the evolution of this structure to the T-shaped case can expect different outputs to the 
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FIG. 23: EM field for the symmetry breaking solution uia/2Tic — 0.3505. Light is incident into the 1 waveguide and scatters to 
the left and right waveguides 2 and 3. P = 0.48W/a. 




FIG. 24: The time dependence of the amplitudes of the light amplitudes \ Ai\, | A2 1 in the cavities (solid and dashed respectively) 
which follow the impulses of the input amplitude Ein (gray). We take the cavities oscillate in non symmetric way: Ai = 0, A2 = 
1. 



right and to the left. Indeed, in spite of the small difference of the defect intensities presented in Fig. [22ja) the 
transmissions and demonstrate vast difference including the case of almost perfect blocking of the transmission 
to the left as shown in Fig. ISSl b) for wa/27rc = 0.3505. Figure [23] shows the anti-bonding BSC is mixed to the 
transport over the input waveguide to give rise to the symmetry breaking solution. 

These results are extremely important for the switching of the output power from the left waveguide to the right 
one. In order to switch the system from one asymmetric state to the other we following Refs. pTl . Iisj apply pulses of 
the input power injected into the waveguide 1. The direct numerical solution of the temporal CMT equation 

iAi = (wi - i"f)Ai + ?:^cr2~e*'^ 

iA2 = (tJ2 - il)A2 + iy/jas-e"^ (63) 

with 5*1+ (t) — Ein{t)e~'^'^* is shown in Fig. [24] which demonstrates the switching effect. The stepwise time behavior of 
amplitude Ein{t) is shown by gray line. One can see that after the first impulse of the input amplitude the oscillations 
of the cavity amplitude relax onto the stable stationary solutions with broken symmetry. Moreover after each next 
impulse the state of the system transmits from one asymmetric state to the other as was observed by Maes et at [ll[ . 



V. DIPOLE MODES OF THE SINGLE NONLINEAR DEFECT COUPLED WITH WAVEGUIDE 



If to present the defect by the single monopole mode i?s(x) we obtain from Eqs. (flT)) 

[lO - LOs - \\As\^ - iVs\As = iy/vlEin, 



(64) 



22 



t = Ei. 



(65) 



where A = —jax'''^H'^s)Es{{^d)'^- That model has attracted interest over the past two decades because of analytical 
treatment and its generality for description of bistability phenomena [l^ |37| - |43| . However the monopole eigen- 
function presents only the trivial identical symmetry transformation. Respectively there is no room for the breaking 
of symmetry. 

In this respect the system becomes nontrivial if we include two eigen dipole modes Ei (x) and E2 (x) of the defect 
rod with the eigen- frequencies uii and W2- The modes are shown in Fig. [25l Respectively, the electric field at the thin 
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FIG. 25: Two defect dipole eigen-modes with the eigen frequencies uji = 0.35 787ra/c and aj2 = 0.36167ra/c in the two-dimensional 
square lattice PhC consisted of the GaAs dielectric rods with radius 0.18a and dielectric constant e = 11.56 where a — 0.5^m 
is the lattice unit. These rods are shown by black open circles. The defect shown by open gray circle has the same radius 0.18a 
and eo = 30. Its center is positioned at Xd = 0,yd = 0.3a. 

defects can be expanded over these dipole modes E{yid) = AiEi{xd) + A2E2{xd) only. Substituting that expansion 
into Eq. (fT6|) we obtain that Eq. (fT7|) will take the following form 



[w - wi - Aii|Ai|2 - A12IA2I2 -I- tT,]Ai - 2X,2ReiAlA2)A2 = iV^{St - S+), 
-2Xi2ReiAlA2 + [lu-lji- A22|^2p - Ai2|Ai|2 + zr2]A2 = iV^{S+ + S+), 



(66) 



where 8^2 fhe amplitudes of light injected simultaneously into both sides of the waveguide. Here with accordance 
to Eq. dill) we have 



16e3/2 



Elix)Elix)d^x. 



(67) 



Moreover we have taken into account the symmetry relations for the coupling constants of the dipole modes with the 
waveguide solutions [1^. As seen from Fig. [25] the hrst dipole mode has the coupling with the left and the right 
ingoing waves opposite signs while the second dipole mode has the same coupling with these waves. Therefore the 
coupling matrix in Eq. (|18p equals 



W ■ 



(68) 



Moreover Eq. (j66p must be complemented by equation for the transmission amplitude 



tLB. = -iVTiAi +i 



tnL = So - i\/TiAi 



T2A2, 

f^A2. 



(69) 



Fig. [26] shows the self-consistent solutions of Eq. (|66)) after substitution of the following model parameters wi = 
0,W2 = 0.01, All = A22 = 0.1, A12 = 0.05, Fi = 0.1, F2 = 0.03. One can see that for the symmetry preserving branch 



23 



7- 
6 

OJ 

< 5 

-;- 4 
< 

3 
2 

1 ^ 
0- 



(a) 






1 

0.8 

1- 






(b) 








'~"o.4 














0.2 























-0.4 -0.2 


0.2 0.4 
o)a/27ic 


0.6 0.8 


0.15 


0.2 0.25 
coa/27ic 


0.3 



FIG. 26: Frequency behavior of (a) intensities of dipole modes and (b) transmissions to the left Ti, = and to the right 

Tr = liiflP for for light injection S^r^ = 0.025 onto both sides of the straight forward waveguide. Blue lines show the intensity 
of even dipole mode, while red lines show the intensity of odd dipole mode, excitement The model parameters are listed in the 
text. 




coa/27tc coa/2jic 

FIG. 27: Frequency behavior of (a) intensities of dipole modes and (b) transmissions to the left and to the right in real 2D PhC 
for 51,2 = 0.015. Dash blue line shows the symmetry preserving solution, solid red line shows the symmetry breaking solution, 
and grey thick line shows the phase parity breaking solution. 



the first even dipole mode is excited with frequency behavior typical for the single nonlinear mode described by Eq. 
(|36p while the second odd dipole mode is not excited because of symmetry. Respectively, the transmission probability 
to both sides of the waveguide equals 1/2 because of normalization condition for the intensity of input light and 
unitarity of the S-matrix \S\\^ + |S'2p + I^laP + l^flLp = 1- However there is also the symmetry breaking branch in 
some narrow frequency domain. Respectively, that branch gives rise to an asymmetry in the light outputs. Moreover 
there is a frequency at which the right (left) output is blocked almost perfectly. 

Also we consider the phenomenon of the symmetry breaking based on the parameters calculated in the 2D square 
lattice PhC with parameters given earlier. The results for the eigen frequencies are collected in Fig. [25] caption where 
the parameters of the PhC are given too. The coupling constants are the following Fi = 0.00075, F2 = 0.00025 in terms 
of I'Kcj a. In order to enhance the coupling constants we substituted two additional linear rods nearby the nonlinear 
defect rod as shown in Fig. [25] Let us evaluate the nonlinearity constants Xmn- With use Ed = eo + 2.yeori2/o, ".2 ~ 

we obtain from 

no ^ 

\nn - -^V^n2h / El^{ic)El{^)cP^. (70) 

We take the linear and nonlinear refractive indexes of the defect rods are, respectively, tlq — = V3, — 

1.8W/a,n2 = 2 X lO^^^cm^/W [H, HI]. Moreover substituting the dipole modes into integrals in Eq. (jTU)) we obtain 
Qii = 0.0123, Q22 = 0.0114, (5i2 = 0.0037. Results of solution of the self-consistency equations for two complex 
amplitudes Ai and A2 are presented in Fig. [^Tfa). Correspondingly from Eq. we obtain the transmission 
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coefficients shown in Fig. [STT b). 

VI. SUMMARY AND CONCLUSIONS 

As by direct solution of tlie Maxwell equations as well as by solution of the CMT equations we have demonstrated 
the symmetry breaking in the system of two nonlinear defects coupled with waveguide through which a light is injected. 
The defects are aligned symmetrically relative to the waveguide as shown in Fig. [1] so there is an inversion symmetry 
relative to the waveguide axis. The thin dielectric rods made from Kerr media are these defects which presented 
by the eigen monopole mode whose eigen frequency belongs to the propagation band of the waveguide. We assume 
that other eigen modes are beyond the band, and therefore have no coupling with the injected light with accuracy of 
evanescent modes. 

That simplest system is remarkable in that it reveals as nonlinearity gives rise to the breaking of symmetry. Indeed, 
let us take temporarily the defects are linear. If the coupled defects were isolated, it would have only two eigen- 
modes, bonding (even) and anti-bonding (odd) with corresponding eigen frequencies (j2ip . For light propagating over 
the waveguide we take that the solution is symmetrical relative to the inversion. Then the light can excite only the 
bonding mode Ag to give rise typical resonance dip at the symmetrical eigen frequency ujg while the anti-bonding 
mode Aa would remain hidden for the propagating light. Obviously, that mode is the simplest case of the bound state 
in continuum .16„ ,67., ,68.J . That state can be superposed to the scattering function with coefficient determined by a 
way to excite the anti-bonding mode 47i] . However the principle of linear superposition is not correct in the nonlinear 
case. Therefore the nonlinearity gives rise to interaction of the scattering state with the anti-bonding mode, i.e., to 
the interaction of the symmetric propagating light with the anti-bonding mode. That obviously breaks the inversion 
symmetry as explicitly shown in Figs. [T^Tb). [TSTc). EOTb). and l^Uf c'). Obviously, the mixing gives total state which is 
nor symmetrical neither anti symmetrical, breaking the mirror symmetry. Further we have shown that the symmetry 
can be broken not only because of different light intensities but also because of different phases of light oscillations at 
the cavities to provoke the Josephson like current between cavities fv^ . For the light transmission in two-dimensional 
PhC the Foyinting power current is an analog of the Josephson currrent, which is shown in Fig. [TW c). It is clear that 
for the light transmission in PhC waveguide there is the power current over the waveguide. However if the waveguide 
is coupled with nonlinear defects the current pattern might be rather complicated even in the symmetry preserving 
scenario as shown in Fig. [HJa) . One can see that laminar current flow over the waveguide induces two vortices around 
the defects which obey the mirror symmetry. As the symmetry has broken the vortex in the waveguide appears which 
is well separated from the defect vortices. Although the mirror symmetry relative to circulation of currents in defect 
vortices is broken, however, the symmetry is remained relative to absolute value of current flows. At last, for the 
phase parity breaking solution still there is symmetry in absolute value of current, but all vortices are exchanging by 
current flows as seen from Fig. [TW c). 

The T-shaped waveguide coupled with two symmetrically positioned nonlinear defects as shown in Fig. [15] can 
be considered as a combination of the previous system and the Fabry-Perot interferometer (FPI) consisted of two 
nonlinear off-channel cavities aligned along the straightforward waveguide considered in Refs. [ll|, [l^, As was 
shown in Ref. [l^l there is a discrete set of the a self-induced bound states in continuum (BSC) which are standing 
waves between off-channel cavities which are to be anti-symmetric in order to elucidate an escape to the input 
waveguide 1. As dependent on position of the nonlinear defects shown consequently in Fig. [T5]the system goes from 
the FPI (case a) to the system considered in section II (case c). In the case (a) the nonlinearity couples the FPI BSC 
shown in Fig. [TH]with the incident wave which is symmetrical relative to the inversion left to right. As the result 
the inversion symmetry breaks as shown for real PhC in Figs. [2TJb) and [2TJc) to give rise to strong asymmetry of 
light outputs to the left and to the right as shown in Fig. mfb). Similarly, the nonlinearity mixes the anti-bonding 
hidden state with the symmetric incident wave as shown in Fig. [53] to give rise different outputs too. At last, for the 
nonlinear defects are coupled as with input waveguide as well as with output waveguides the frequency behavior of 
the transmissions is complicated as shown in Fig. [22lb'). In Ref. [H^ we demonstrate as these phenomenona can be 
explored for all-optical switching of light transmission from the left output waveguide to the right one by application 
input pulses. 

In section V we show that breaking of symmetry might occur even for the single nonlinear defect positioned in the 
straight forward waveguide provided that the defect is presented by two dipole modes. That model embraces two 
degrees of freedom of the defect while in the former models each defect was presented by the only degree of freedom. 
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